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Abstract 



PACS: 12. 38. Aw, 12.39.Ki, 12.39.Pn 



Two-dimensional QCD is studied from the point of view of existence of multiple 
chirally noninvariant solutions to the mass-gap equation. The ground-state solution is 
reproduced and an infinite set of replica solutions is discovered for this equation using 
the WKB quantisation procedure. 
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In this paper, we revisit the two-dimensional model for QCD [UI2] in the axial (Coulomb) 
gauge [Him an d address the question raised in a set of recent publications [SUHllZj as if replica 
solutions to the mass-gap equation exist defining metastable chirally noninvariant phases of 
the theory. For the four-dimensional potential model for QCD (HI E] such solutions were 
discovered, for the oscillator- type potential, in [Hj, for the phenomenological potential with 
the Coulomb interaction included, in [S], and, finally, for an arbitrary power-like potential, in 
A similar conclusion was also made in a different approach [TUj. In this paper, we apply 
the analytical and numerical techniques, developed in the work [7j for the aforementioned 
four-dimensional model for QCD, to a simpler theory — the 't Hooft model for QCD in 
two dimensions, and demonstrate how replicas appear in this case. Thus we reproduce the 
numerical solution for the BCS vacuum of the theory found in [IT] and build explicitly, as 
an example, the first "excited" solution — the replica — for the mass-gap equation. We 
argue that the number of such solutions is infinite and evaluate their behaviour at the origin, 
considering the mass-gap equation as a Shrodinger-like equation and using the quasiclassical 
quantisation procedure for the latter. 

Let us give necessary details of the two-dimensional model for QCD suggested by 't Hooft 
in 1974 £Q. The theory is given by the Lagrangian 

L(x) = -±F; u (x)F; u (x) + q(x)(i^d, - gA$* 7lt - m)q(x), (1) 
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and the large- Nc limit is understood, such that the effective coupling constant remains finite, 



9 2 N C 

7 = — — ► const. (2) 

It was demonstrated in ^2] that, in such a limit, the well-known no-go Coleman theorem 
[TB*] can be bypassed, and the spontaneous breaking of the chiral symmetry (SBCS) may 
happen, provided the limit Nq — > oo is taken prior to the chiral limit, m — > [T3] . 

We gauge the theory by the condition 

A 1 (x o ,x) = 0, (3) 

so that the only nontrivial correlator of the gluonic fields, the gluon propagator, takes the 
form: 

D$(x -y ,x-y) = - % -5 ah \x - y\6(x - m), (4) 

yielding a linear confinement for the inter-quark interaction mediated by the two-dimensional 
gluon. We use the following convention for 7 matrices: 

70 = [3 = 03, 71 = ia 2 , 75 = ol = 7071 = a u (5) 

cr's being the standard 2x2 Pauli matrices. It was demonstrated in [3] that, in the axial 
gauge Q, the infrared divergences in the theory can be self-consistently regularised by 
the principal value prescription, which is understood in all integrals below. If not stated 
otherwise, we consider the chiral limit, m = 0. 

Following the standard technique, we build the Hamiltonian of the theory, 

f 9 o 2 f 

H = ~ j dxq+ ( tj x ' )75 a^ <? ^' ^ _ ~4 J dxdy q ^ x ) ta ^> x ) \ x ~ y\ « (*> s/)*°0(*> y)> ( 6 ) 

and arrange the normal ordering of the latter in terms of the dressed quark field, 

q ai (t,x) = j ^e ikx [b a (k,t) Ul (k) + di(-k,t)v t (-k)}, (7) 

{b a (t,p)bl(t,q)} = {d a (t,-p)dl(t,-q)} = 2ir5( P -q)5 ct p, (8) 
1\ , > „, > /o 



u(p) = T{p) (\ , v(-p) = T(p) M , Tip) 



cxp 



(9) 



The Bogoliubov (chiral) angle 6 P (shorthand notation for 9{p)) is odd and 9{p — > ±00) — > 
±7r/2. The normally ordered Hamiltonian © splits into the vacuum energy, the quadratic, 
and quartic parts in terms of the quark creation/annihilation operators, 

H = LN C £ V + : H 2 : + : H 4 :, (10) 

where L is the one-dimensional volume; the vacuum energy density reads: 

£v = I 2% Tr {™ A -(P) + H Jp-^ A +^ A -^} ' = TiP) l -^T\p). (11) 
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It is convenient to define the excess of the vacuum energy density over the free case 9f vee (p) 
fsign(», 



AS V [9) = S v [9] - S v [9 iTCC ] = -J ^( P sin9 p -\p\)-l-j cos[9 p - 9 k \. (12) 



dpdk 
(p - k:f 



The actual form of the chiral angle is such that the quadratic part of the Hamiltonian 
(jlOjl becomes diagonal in terms of the quark creation/annihilation operators, or alternatively, 
the vacuum energy takes its minimal value: 



which yields the mass-gap equation for the chiral angle 9 P [3], 



pcosi 



6k] 



(13) 



(14) 



where, as mentioned above, the integral is defined by means of the principal value prescrip- 
tion. For a given solution of Eq. (JT3J), the dressed quark dispersive law can be found from 
the system of equations [3]: 



_ 7 f dk 



dk 



cos 9i 



E p sin 9 P 



p + 



and reads 
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2 7 (p - k) 2 

dk 



(15) 



sin 9k 



E p = psm9 p + --J- - — - cos[9 p - 9 k 



The mass-gap equation (|14j) is the main object of our studies. Following the method 
suggested in [7j, we rewrite (|14j) in the form of a Shrodinger-like equation for the wave 
function ip p = cos 6^, in coordinate space, with the help of the relation ()16|) : 



- 7T7 

E„ + — \x 



Jp ' 2 

where the following Fourier transformation was used: 



0. 



(17) 



I 



dk 



-7T 



dx\x\ip x e ipx . 



J Jp^kJ^ k 

For further analysis it is convenient to introduce the effective mass of the quark m„, 



sin 6L 



V 



COS 9 r , 



m r . 



p 2 + m 2 ' Jp 2 + m 2 



(19) 



and to rewrite the dispersive law E p in terms of the latter 1 : 



E„ 



_Xf 



dk 



kjp 2 + m 2 



P +mp 2l(p-k) 2 [ pJk^ 2 



(20) 



1 Notice that Eqs. I|17|) and l|2(J|) follow immediately from the system (|15H . form the first and the second 
equations, respectively, and using the fact that, by definition of the principal value prescription, f (p^k)' 1 = ^- 



so that 

E P « \p\ - (21) 

where the approximation m p <C p used in the last formula is valid for momenta p's larger 
than some scale po generated by the function m p , po ~ m(p = 0). Such a simplification is 
justified by the form of the function m p which takes its maximal value at p = and then 
decreases fast for finite momenta. It is clear, however, that in the approximate formula (|2*T|) 
we suppressed the positive contribution coming from the first term on the r.h.s. of Eq. (J2UJ) 
and enhanced the negative contribution coming from the second term on the r.h.s. As a 
result, the approximate mass-gap equation takes the form: 



% = (22) 



where we passed over to the momentum space, considering the quark dispersive law E p as an 
effective potential. We are interested in the negative eigenvalues e. Besides that, according 
to the definition of the wave function, ip p = cos# p , only even eigenstates of the equation (J22I) 
should be selected. We use the quasiclassical quantisation procedure to find the spectrum of 

Eq. U22D, 

vr [in + -) , (23) 



'm ax 



7T7 JO 



P + £ 

p 
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and see that the integral on the l.h.s. diverges at p = since the approximation made above, 
m p <C p, fails for small momenta, where m p tends to a nonzero limit m(p — 0) = fi. Therefore 
we introduce fi into the integral (J23~j) . as a cut-off, for example, through the substitution 
\p\ — > \/p 2 + A* 2 in the quark dispersive law. The resulting potential V(p) = E p appears 
sufficiently attractive to possess bound states. Notice, however, that the cut-off /i itself is not 
a constant, but it decreases fast for each next solution. Moreover, the solutions of the exact 
mass-gap equation (fTf]) correspond to e n = 0, for any n, so that one can reverse the arguments 
and use this condition to evaluate the behaviour of [i n for an arbitrary n. The regularised 
quasiclassical integral, with e put to zero, is easily done, /wkb ~ — - ^{^n/ ^/l) + const, 
yielding 



TT 2 



Mn = \/7expl — -n + 5\, n = 0,1,2,..., (24) 



where S is an n-independent finite correction to the leading regime. From Eq. (|24j) one can 
easily estimate the behaviour, at the origin, of the solution 9 P to the mass-gap equation (fTlj) . 
Indeed, 

9 P ^ sm6 p ^ — = const— e^ 2/2)n , (25) 



where, in order to find the constant on the r.h.s., one needs either to find at least one 
solution to the mass-gap equation numerically, or to estimate it using a more accurate 
regularisation procedure for the dispersive law and evaluation of the quasiclassical integral. 
From the numerical solution for the ground-state vacuum found in ^TJ |3] this constant 
equals, approximately, 1.2. The accuracy of the approximate formula (|23j) is expected to be 
better that 10% even for the lowest states, as it follows from the consideration in ref. [7j. 

In Fig. 1, as examples, we give the profiles of two solutions to the mass-gap equation 
(H3J), found numerically, corresponding to the BCS vacuum (6* ) and to the first replica (9i). 
The behaviour of these two solutions is in agreement with the general formula (J2SJ). 
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Figure 1: The two lowest solutions to the mass-gap equation (jUj) corresponding to the BCS 
vacuum (n = 0) and to the first replica (n = 1). Since the chiral angle is an odd function, 
we plot it only for positive momenta. Momentum p is given in the units of y/j. 

From the formula (|25J) one can make several conclusions: i) the solution (|25j) is purely 
nonperturbative since it involves y/j, whereas only integer powers of 7 appear in any per- 
turbative expansion; ii) the mass-gap equation of the 't Hooft model contains an infinite 
number of solution — replicas — which appear due to a very peculiar sharp behaviour 
of the quark self-energy at small momenta, where the second, negative, term on the r.h.s. 
in Eqs. (120)1 . (|2*T|) allows fast oscillations of the chiral angle in the vicinity of the origin. 
Each replica defines a phase of the theory with spontaneously broken chiral symmetry; iii) 
the parameter defining SBCS, for example fj, n , decreases very fast for each next solution, 
/i n+ i//i n ~ exp (— 7r 2 /2) ~ 10 -2 , which makes excited solutions extremely hard to find nu- 
merically; iv) the series of replicas, for n — > 00, converges to the trivial chirally symmetric 
solution trivial (p) = f sign(p), with a discontinuity at p — 0, which can be formally viewed 
as the limiting case n = 00 of the formula (|2"H|) : v) the vacuum energy density of repli- 
cas increase for each next solution converging to the trivial vacuum energy. In contrast to 
four-dimensional QCD, where the trivial vacuum corresponds to zero vacuum energy, in two 
dimensions the trivial vacuum energy density is infinite (see the discussion in 

As discussed in the papers [2J E] , the existence of replicas is rather a general feature of 
confining interactions, at least as far as confining potentials are concerned. In ref. [7] a similar 
consideration was made for a generalised power-like potential in four dimensions and the 
same conclusion was made concerning an infinite number of solutions to the corresponding 
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mass-gap equation, as a consequence of a very peculiar behaviour of the quark self-energy 
in the infrared domain. Still it is instructive to give an example of the model where there 
exists a unique solution to the mass-gap equation. Using the same technique as in ref. jjj, 
one can easily generalise the formula (J25|) for the case of an arbitrary power-like potential in 
two-dimensions, V(x) 



K« +1 \x\ 



6„ oc 



P 



, p ^ ■ exp 



vrn | r Q) r ( V' 



l/a 



r(i±a) 



(26) 



which reproduces Eq. (|25j) for a = 1; r(^) being the Euler Gamma-function. One can see 
that, for a — > 2, that is, for the harmonic oscillator-type potential, the exponential factor on 
the r.h.s. of Eq. (|2*Bj) becomes infinite, so that only the case of n = — the BCS vacuum 
- survives. In order to illustrate the situation better let us apply the general method 
used in this paper to the harmonic oscillator-type potential directly. First, we derive the 
formula for the quark self-energy for such an interaction. Using Eq. (|16|) . with the change 
l/ip ~ k) 2 — > ttKq5"(p — k) (see j2] for the details), one arrives at 



E p = p sin 6 P -^K*[6. 



T 
pi 



71" 



\P\ - y*', ! [*(p)P 



where, in the last formula, we used the well-known representation for the delta-function, 
S(p) = lim^o ~ ^2^~2 • As a result, the mass-gap equation, similar to (J22J), becomes a second- 
order Schrodinger-like differential equation, 



-iKl " 



dp' 



(28) 



but the corresponding Bohr-Sommerfeld integral, for vanishing eigenenergy, monotonously 
increase from zero to its maximal value of n, reached for jj, — > 0. This makes a crucial 
difference between this theory and the 't Hooft model, with linear confinement, where the 
Bohr-Sommerfeld integral is not bounded from above and grows logarithmically with the 
decrease of the regulator fi. Therefore, for the two-dimensional linear confinement, the WKB 
quantisation condition can be satisfied for an arbitrary n, provided fi n behaves properly (see 
Eq. (125)1 above), whereas for the harmonic oscillator-type potential this condition can be 
fulfilled for the only value n = 0. Thus we repeat the conclusion deduced from the general 
formula ([26)1 . Indeed, the mass-gap equation ([28)1 . which can be rewritten in the form 



p cos Or, 



(29) 



is known to have a unique nontrivial solution which can be found either numerically 0, or 
it can be built approximately using the method of separation of large- and small-momentum 
regions where the simplified equations are solved analytically with the consequent matching 
of the solutions found in different intervals (see jH] for the four- dimensional case). The result 
reads: 

' 2 W^) 1/3 Mj sign(p), (30) 



6n ~ — 

p 2 



1 - 



32/3r Q 



with Ai(z) and T(z) being the Airy and Gamma-function, respectively. The formal reason 
for this behaviour of the theory can be easily deduced, as before, from the form of the quark 
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self-energy (j27j). where the negative term has a delta-functional form and contributes only at 
p = 0, whereas, in order to support replicas, it must give a sizable effect at finite momenta. 

In four-dimensional QCD, the mass-gap equation for the harmonic oscillator-type poten- 
tial is also differential, but it is known to support an infinite tower of replicas [S], similarly 
to other power-like potentials in four dimensions [7j. Extra singular terms in the quark self- 
energy appear, in this case, due to a reacher structure of the interaction in four dimensions, 
namely, due to presence of the angular variables. The formula similar to (f2*U|) is derived in 
ref. |7] and it contains an extra factor which cancels the singularity at a = 2 and makes the 
corresponding limit smooth — namely, this formula coincides with Eq. (126)) with the only 
change T (pp) — ► T (pp) [I], so that the general expression can be written with T (pp), 
D being the number of dimensions. The case a = D is degenerate, that is, a unique solution 
to the mass-gap equation exists, but it can be reached only in two dimensions since, as 
argued in for a > 2, the mass- gap equation is divergent. 

In contrast to four- dimensional QCD, where the limit Nc —>■ oo can be relaxed and the 
effect of the QCD string breaking is expected to reduce the number of replicas making it 
finite f7J, in two dimensions, as mentioned above, the requirement Nc = oo is necessary to 
bypass the Coleman theorem and allow SBCS. 

Let us discuss a possible interpretation of replicas and their effect on two-dimensional 
observables. Notice that a whole tower of single-quark states can be built over each vacuum 
by a repeated application of the corresponding creation operators to the vacuum. Any such 
set provides a complete basis sufficient to saturate any multiquark state in the theory. An 
alternative, and in many cases more convenient, mathematical language for such sets is 
given by the well-known coherent states. Indeed, any nontrivial vacuum of the theory \n), 
n = 0,1,2..., can be obtained from the trivial vacuum |0)o, corresponding to the chiral 
angle 9 p m = |sign(p), by means of the unitary operator S n [To] : 

|n) = S n \0) , S n = exp (q{ -Q n ), QI = \eJ ~ ^ riv )^)4(-p), (31) 

or in the non-exponential form: 

(A9 n A9 n \ 

cos^ + Ct a sin^J|0) , A9; = 9 n p -9^\ C\ a = &t(p)4(-p). (32) 

It is easy to derive the vacuum orthogonality condition now, which reads 



On Qm 

(m\n) = o(0|^S n |0) = II cos ^- ~ = ex P 



LNc / In cos — ?— 




0, (33) 



where L is the one-dimensional volume. Thus any two vacua become orthogonal in the 
infinite volume. 

Let us define coherent-like states \<p,P) 0: 

\<p, (3) = J] (cos ^ + ClM sin ^) |0) = II \<P, P)p*, ClM = e^ sign(p) Ct, (34) 

p,a \ ^ & / p t a 

which form an overcomplete set in the Fock space of the colourless quark-antiquark pairs, for 
any given momentum p spanned by the two orthogonal states, {|0) , Cp a [/^p]|0)o}- Therefore 
one can resolve the unitary operator in this space in terms of the states (|34|) : 

dip pa / d/3 pa \tp,/3) pa J pa (tp,a)pa((p,{3\. (35) 

... -IT JO 



p,a 
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where the choice of the Jacobean J is not unique, the simplest form being just 2ir 2 . The BCS 
vacuum state and all replicas belong to the set ([34)1 and correspond to a specific choice of the 
parameters ip and (3, \n) = \A9 n , |). Therefore the coherent-like states (J34j) define an over- 
complete basis, alternative to the complete set of multiquark states, each containing a fixed 
number of particles. Any vacuum replica can be chosen as the ground state to produce such 
an overcomplete basis, and all other replicas look like the coherent-like states in this basis. 
Notice that such a picture is maintained when the Hamiltonian of the theory is truncated 
and the : H4 : part in the expression (JIUj) is neglected. Meanwhile as soon as the quartic 
part : H4 : is included, the appropriate degrees of freedom of the theory, which minimise 
the vacuum energy and diagonalise the full Hamiltonian, are given by the quark-antiquark 
mesons. A second, generalised, Bogoliubov transformation is a possible way to arrive at 
the fully diagonalised Hamiltonian and at the bound-state equation, which plays the role 
similar to the mass-gap equation for the quadratic part of the Hamiltonian . The hadronic 
(mesonic) spectrum of the theory constructed over any replica vacuum contains tachyons, 
only the true BCS vacuum being free of the latter. Therefore the interpretation of replicas as 
stable vacua does not hold anymore, whereas, providing local minima of the vacuum energy, 
they may realise as metastable states in hadronic reactions j^j. In what follows we confine 
our consideration of the trivial vacuum |0)o, so that for the sake of simplicity, we use the 
notation A9™ for the difference between the chiral angle for the nth replica |n) and the BCS 
vacuum |0). 

In the paper jH] a field theory approach to replicas in four dimensions is suggested which 
allows one to incorporate replicas, as local in space-time objects, into quark models. The 
given approach is easily adapted for the two-dimensional 't Hooft model, so that the local 
operator creating the nth replica reads 0: 



s n = Tex P 



d 2 xf {n \x )ip(x ,x) (-^75-7-) ip(x ,x) 



dx / 

where the function f^ n '{xo) is even with the Fourier transform being such that 



(36) 



A9 n 

f<»\2E p ) = (37) 

On substituting the representation (JJJ) for the quark field to Eq. (|36jh and with the help of 
the definition (|3*7j) . one easily recovers the form (J3l"|) for the operator S n . 

From the operator (|3T)|) one can easily derive the diagrammatic technique involving the 
new quark-quark-replica vertex, and the orthogonality condition (|3*3*|) can be easily repro- 
duced then, via summation of an infinite set of fish-like diagrams, as a simple selfconsistency 
test of the approach. Performing a systematic expansion of all expressions in terms of the 
small difference A9™ and truncating the series at the second power of the latter, one can 
arrive at a simplified diagrammatic technique with an effective propagator of the replica: 

G n (p ,p) = u n f^( Po )f( n \- Po )(2ir)6(p), (38) 

where u n = N n /T, T being the full time of the process, N n giving the number of times 
the nth replica was excited during the time T. This propagator is attached to the quark 
lines through the two vertices An overall sum over n is understood for the amplitude 
of the process in order to take into account the contribution of all replicas. In hadronic 
processes, the diagrams with replicas excitation provide corrections to standard amplitudes. 
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This work is in progress now and will be reported elsewhere. The interested reader can find 
further details of the formalism in the paper 0. To conclude, let us mention that with the 
replica solutions to the mass-gap equation found in this paper, the 't Hooft model confirms 
its status of a marvelous test-ground for the four-dimensional QCD where many interesting 
phenomena can be observed and various approaches can be checked and tuned. 

One of the authors (A.V.N.) would like to acknowledge the financial support of IN- 
TAS grants OPEN 2000-110 and YSF 2002-49, as well as the grant NS-1774.2003.2, and 
the Federal Programme of the Russian Ministry of Industry, Science and Technology No 
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